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Abstract 



We show that a natural realization of the thermostatistics of g-bosons can 
be built on the formalism of q-calculus and that the entire structure of ther- 
modynamics is preserved if we use an appropriate Jackson derivative in place 
of the ordinary thermodynamics derivative. This framework allows us to ob- 
tain a generalized g-boson entropy which depends on the g-basic number. 
We study the ideal g-boson gas in the thermodynamic limit which is shown 
to exhibit Bose-Einstein condensation with higher critical temperature and 
discontinuous specific heat. 

PACS numbers: 05.30.-d; 05.20.-y; 05.70.-a 
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I. INTRODUCTION 



The spin-statistics theorem represents one of the fundamental principles of physics and 
establishes a strict connection between quantum mechanics of many-body systems and quan- 
tum statistical mechanics. The complete symmetrization or antisymmetrization of the many 
body wave function (or the commutation-anticommutation relations in the language of sec- 
ond quantization) reflects the contrasting nature of bosons and fermions. Such quantum 
many body statistical behavior affects the number of possible states of the system cor- 
responding to the set of occupation numbers and consequently the collective statistical 
mechanics description. 

The power of the statistical mechanics lies not only in the derivation of the general laws 
of thermodynamics but also in determining the meaning of all the thermodynamic functions 
in terms of the microscopic interparticle interaction and in providing a collective description 
of the equilibrium many body system by means of the macroscopic variables such as pressure 
and internal energy. 

In the recent past there has been increasing emphasis in quantum statistics different 
from the standard bosons and fermions. Since the pioneering work of Gentile and Green 
IIHIJ], there have been many extensions beyond the standard statistics (such as parastatis- 
tics, fractional statistics, quon statistics, anyon statistics and quantum groups) which have 
become topics of great interest because of the wide range of applications envisaged, from 
cosmic strings and black holes to the fractional quantum Hall effect and anyonic physics in 
condensed matter [||]. 

In the literature there are two principal methods of introducing an intermediate statistical 
behavior. The first is to deform the quantum algebra of the commutation-anticommutation 
relations thus deforming the exchange factor between permuted particles. The second 
method is based on modifying the number of ways of assigning particles to a collection 
of states and thus the statistical weight of the many-body system. The two methods are 
related but a full connection between the quantum mechanics approach and the statistical 
mechanics approach is possible only with the simultaneous knowledge of both. 

One interesting realization of the first approach is the study of exactly solvable statistical 
systems which has led to a new algebra, the g-deformed algebra of creation and annihilation 
operators, usually called g-bosons (g-fermions) or g-oscillators and related to the general 
theory of quantum groups Many recent investigations in the theory of g-bosons have 

provided much insight into both the mathematical development and the g-deformed ther- 
modynamics |6|-|l4}j . However, we believe that a fully consistent formulation connecting the 
statistical mechanics and the thermodynamics (i.e., thermostatistics) of g-bosons has been 
lacking. In particular it is desirable to derive an explicit expression for the entropy of the 
g-bosons, which plays a central role in the thermostatistics of the system and in the informa- 
tion theory. It is important to show that the full structure of thermodynamics of g-bosons 
is preserved and the closed loop of thermodynamic relations is satisfied. This is a nontrivial 
task because there is no a priori reason that the thermodynamic relations be automatically 
preserved for the g-deformed structures. 

A remarkable example is the Tsallis nonextensive statistics [(L5[], based on a generaliza- 
tion of the Boltzmann-Gibbs entropy, where the thermodynamic functions such as entropy 
and internal energy, are deformed but the whole structure of thermodynamics is preserved. 
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Although Tsallis nonextensive thermodynamics is inspired by the (multi)fractal property 
of a system and does not embody quantum group theory, many papers are devoted to the 



formal analogies between (/-oscillators and nonextensive statistics [|I6|-|T9|| . The reason for 



this connection has to do with the common language of the two deformed theories which is 
the g-calculus. 



The g-calculus was introduced at the beginning of this century by Jackson |20| in the 
study of the basic hypergeometric function and it plays a central role in the representation of 
the quantum groups |^TJ . In fact it has been shown that it is possible to obtain a "coordinate" 
realization of the Fock space of the g-oscillators by using the deformed Jackson derivative 
(JD) [^,[23|]. Moreover we observe that it has recently been shown that the JD can be 



identified with the generators of fractal and multifractal sets with discrete symmetries p4 



Since the thermodynamic functions of nonextensive statistics are deformed by using the 
framework of g-calculus, we expect g-calculus to play an important role also in g-boson 
thermostatistics. 

It is the purpose of this paper to show that a fully consistent thermostatistics of q- 
boson gas can be obtained by using an appropriate Jackson derivative rule in the standard 
thermodynamics relations. In this framework, the whole structure of thermodynamics is 
preserved and this enables us to derive all the thermodynamic quantities such as the entropy, 
internal energy and the distribution function in the g-deformed theory Special attention is 
paid to the study of the ideal g-boson gas and the phenomenon of g-boson condensation. 

This paper is organized as follows. In Sec. II we review the g-boson algebra and outline 
the modification of the standard boson theory brought about by the g-calculus. In Sec. Ill 
we determine the distribution function of the g-boson gas by utilizing the standard definition 
of the thermal average of an operator. In Sec. IV we introduce a consistent prescription for 
the use of the Jackson derivative in the thermodynamic relations. This allows us to obtain in 
Sec. V the generalized entropy for g-bosons and to derive this from the deformed statistical 
weight. Sec. VI describes the behavior of the ideal g-boson gas and the phenomenon of 
g-boson condensation. We report our conclusions in Sec. VII. 



II. Q-BOSON ALGEBRA AND ITS REALIZATIONS 

We shall briefly review the theory of g-deformed bosons defined by the g-Heisenberg 
algebra of creation and annihilation operators of bosons introduced by Biedenharn and 
McFarlane 0,0], derivable through a map from SU(2) q . The g-boson algebra is determined 
by the following commutation relations for a, and the number operator N, thus (for 
simplicity we omit the particle index) 

[a, a] = [a\ a'] — , aa) — qa)a = 1 , (1) 

[N, a f ] = a f , [N, a] = -a . (2) 

The g-Fock space spanned by the orthornormalized eigenstates |n) is constructed accord- 
ing to 
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|n> = -y=|0), a|0) = 0, (3) 



where the g-basic factorial is defined as 

[n]\ = [n][n-l] • • • [1] (4) 
and the g-basic number [x] is defined in terms of the g-deformation parameter 

H . ^1 . (5) 

For the following discussion it is worth observing that the g-basic number satisfies the 
non-additivity property 

[x + y] = [x] + [y] + (q - 1) [x] [y] . (6) 

In the limit q — > 1, the g-basic number [x] reduces to the ordinary number x and all the 
above relations reduce to the standard boson relations. 
The actions of a, on the Fock state \n) are given by 

aV) = [n+l] 1/2 |n+l> , (7) 
a\n) = [n] 1/2 \n- 1) , (8) 
N\n) = n\n) . (9) 

From the above relations, it follows that a^a = [N], aa^ = [N + 1]. 

We observe that the Fock space of the g-bosons has the same structure as the standard 
bosons but with the replacement n\ — > [n]\ . Moreover the number operator is not a^a but 
can be expressed as the nonlinear functional relation iV = f(a'a) which can be explicitly 
written formally in the closed form 

N = -!— log (l + (g - l)a f a) . (10) 
log g v ' 

The transformation from Fock observables to the configuration space (Bargmann holo- 
morphic representation) may be accomplished by choosing |22,23| 



a) = x , a = d { x q) , (11) 
where is the Jackson derivative (JD) |2(J 

^./MzM, (12) 

which reduces to the ordinary derivative when q goes to unity and therefore, the JD occurs 
naturally in g-deformed structures [Ell]. 
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III. THERMAL AVERAGES AND STATISTICAL DISTRIBUTION FOR 

Q-BOSON GAS 



Several investigators have studied the equilibrium statistical mechanics of the gas of 
non- interacting g-bosons [|BHl~4||. We shall now briefly discuss some of the important results 



from these studies before introducing our formulation of the thermostatistics of g-deformed 
bosons. 

In the grand canonical ensemble, the Hamiltonian of the non-interacting boson gas is 
expected to have the following form 



H = '£(e i -n)N i , (13) 

i 

where the index i is the state label, /i is the chemical potential and is the kinetic energy 
in the state i with the number operator iVj. It should be mentioned that the form of the 
Hamiltonian is not unique in the literature, where some authors introduce the Hamiltonian 
which involves the basic number [iVj]. The advantage of the form in Eq.flTB"!) is that it 
describes clearly the number of particles in the spectrum by an integer number and will 
allow us to generalize the laws of thermodynamics in a simple manner. 
The thermal average of an operator is written in the standard form 

Tr (Oe-P H 



(O) = — ^ , (14) 

where Z is the grand canonical partition function defined as 

Z = Tr(e- pH ) , (15) 



and f3 = 1/T. Henceforward we shall set Boltzmann constant to unity. Let us observe that 
the structure of the density matrix p = e~^ H and the thermal average are undeformed. As a 
consequence, the structure of the partition function is also unchanged. We emphasize that 
this is not a trivial assumption because its validity implicitly amounts to an unmodified 
structure of the Boltzmann- Gibbs entropy, 

S = logW, (16) 

where W stands for the number of states of the system corresponding to the set of occupation 
numbers {rii}. Obviously the number W is modified in the g-deformed case. It may be 
pointed out that in the case of nonextensive g-deformed Tsallis statistics, the structure of 



the entropy is deformed via the logarithm function [15 ]. 

By using the definition in Eq. @ of the g-basic number, the mean value of the occupation 
number ni can be calculated starting from the relation 



i Tr (e-^ajoi) , (17) 



L m 2 

and after applying the cyclic property of the trace and using the g-boson algebra, it is easy 
to show that 
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'"" e -^-M) . (is) 



[m + 1] 



The explicit expression for the mean occupation number can be obtained by using the 
following property of the basic number, 

[m + l] = q N + 1 , (19) 

and hence for q real, 

r^^log * (20) 
logg \z- L eP^ — q J 

where z = is the fugacity. It is easy to see that the above equation reduces to the 
standard Bose-Einstein distribution when q — ► 1. The total number of particles is given by 
N = Ei rii. 



IV. JACKSON DERIVATIVES IN Q-THERMODYNAMICS RELATIONS 

From the definition of the partition function, Eq. flT5|) , and the Hamiltonian, Eq.(^), it 
follows that the logarithm of the partition function has the same structure as that of the 
standard boson 

\ogZ = -£log(l -ze^) . (21) 

i 

This is due to the fact that we have chosen the Hamiltonian to be a linear function of 
the number operator but it is not linear in a^a as seen from Eq. (|T0|) . For this reason, the 
standard thermodynamic relations in the usual form are ruled out. It is verified, for instance, 
that 

N^z^-logZ. (22) 

As the coordinate space representation of the g-boson algebra is realized by the introduc- 
tion of the JD (see Eq.(|Tl|)), we stress that the key point of the g-deformed thermostatistics 
is in the observation that the ordinary thermodynamics derivative with respect to z, must 
be replaced by the JD 

l z ^ , (23) 

where we have defined as the Jackson derivative up to a constant (which goes to unity 
when q — > 1) 

= Izl q(m) . (24) 
logg 

Consequently, the number of particles in the g-deformed theory can be derived from the 
relation 
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N=zV^\ogZ = Y, n i > 

i 



(25) 



where n« is the mean occupation number expressed in Eq.(|20|). 

The usual Leibniz chain rule is ruled out for the JD and therefore derivatives encountered 
in thermodynamics must be modified according to the following prescription. First we 
observe that the JD applies only with respect to the variable in the exponential form such 
as z = or y t = e _/3e \ Therefore for the g-deformed case, any thermodynamic derivative of 
functions which depend on z or must be converted to derivatives in one of these variables 
by using the ordinary chain rule and then applying the JD with respect to the exponential 
variable. For example, the internal energy in the g-deformed case can be written as 



E^f£'log(l-^<). (26) 



In this case we obtain the correct form of the internal energy 

U = '£e i n i , (27) 



where n, is the mean occupation number expressed in Eq.(|20|). 

This prescription is a crucial point of our approach because this allows us to maintain 
the whole structure of thermodynamics and the validity of the Legendre transformations in 
a fully consistent manner. 



V. ENTROPY OF THE Q-BOSON GAS AND THE DEFORMED STATISTICAL 

WEIGHT 



In light of the above discussion, we have the recipe to derive the entropy of the g-bosons 
which leads to 



S 



df 



dKj 
8(5 



V$ log(l - 



K i 



log Z + (3U - (3 fiN 



where Ki = ze~^\ U and N are the modified functions expressed in Eqs. ([ 
Q = —TlogZ is the thermodynamic potential. 

Using Eqs.(|18D-([2"0D, after some manipulations, we obtain the entropy as follows 

S = { ~ n i lo § \ n i\ + ( n i + !) lQ g \ U i + !] _ U i lo S l} ■ 



(28) 

and (|25T) and 



(29) 



The above entropy goes over to the standard boson entropy in the limit q — *■ 1. It has the 
compact form which resembles the entropy of the standard boson but with the appearance 
of the g-basic numbers, [ni\ and [n^ + 1], in the argument of the logarithmic function and in 
the presence of the last term, — n^logg, which follows from non-additivity property of the 
g-basic number. In fact, using Eqs.(^) and (|I9"D can be re-expressed as 
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Hi log q = log ([rii + 1] - [rn]) 



(30) 



The expression for the entropy is very relevant to the statistical information about the 
number of possible states occupied by the g-bosons and gives us the desired connection 
between the deformed quantum algebra and the quantum statistical behavior. It is inter- 
esting to observe that in the classical limit, the entropy does not reduce to the standard 
Boltzmann-Gibbs entropy (S = — J2i n>i logn,), but remains deformed, except in the limit 
q — ► 1. This result is similar to the case of Greenberg's infinite statistics and the quantum 



Boltzmann distribution obtained as a particular case of quon statistics [25]. The meaning 



of this is that the deformation exhibited in the entropy transcends the quantum nature but 



is built into the theory, somewhat similar to the case of nonextensive Tsallis statistics [JT5 
The origin of the connection between the two different deformations (g-deformed quantum 
groups and nonextensive statistics) is beyond the scope of this paper and will be reported 
elsewhere. 

In order to assure consistency, we must now show that the extremization of the entropy 
with fixed internal energy and number of particles leads to the correct g-boson distribution 
function. The extremum condition can be written as 

8 (S - /3U + /3/iiV) = , (31) 

where (3 and (3li plays the role of Lagrange multipliers. 

To perform such extremization in the g-boson case, we assume that the mean occu- 
pational number depends on the energy only as a function of jji = e _/3e % S = S[n(yi)]. 
Following our prescription described in Sec. IV on the use of JD, the above extremization 
condition can be written as 

V$(S-pU + l3tiN)6y i = 0. (32) 



Employing Eqs. (25), (p7|) and (p9j), and carrying out the JD, the extremization condition 
reduces to 

<m) ( log MM±il _ A - n(Vi) ( log If!M±i! _ A + 
V q[n(m)] J \ Hvi)] J 

n( yi )logq-log ^1^1 =0, (33) 
[n{qyi) + 1] 

where 6j = /3(ej — li). 

The algebraic simplification of the above equation is intractable because of the complexity 
of the property of g-basic numbers. However, it is possible to determine the solution of the 
equation by observing that for any function f(x), there exists a functional relationship 

q " \f(qx) + 1] [f(qx)\ q [f(x)\ • 1 ] 

The notation <^==^> used here denotes that one relation implies the other and vice versa. The 
validity of the first relation in the above equation eliminates the last two terms in Eq.(p3|) 
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and the validity of the second relation implies that the quantities in parenthesis in Eq.(|33|) 
are equal, and since n(qyi) ^ n(yi), for q 1, it follows that Eq.(^) is satisfied if 

Hm) + c* . (35) 



Hvi)) 

The above relation is equivalent to Eq. (|1|) which implies the mean occupational number 
Hi of Eq.®. 

As discussed earlier, the entropy provides the information about the statistical weight 
W which will be deformed in the case of g-boson particles. To investigate this deformation 
we begin with the basic relation for the entropy 

S = \ogW q , (36) 

where W q is the deformed statistical weight. Just as the ordinary factorial n\ is replaced by 
the g-basic factorial [n]! in the construction g-Fock space (see Eq.(||)), we assume that this 
substitution also prevails in the expression for the statistical weight and hence we require 

W -T\h±MzJ$ (37) 
where g± denotes the number of subcell levels. The reason for this modification lies in the 



definition of the binomial coefficient in the g-combinatorial calculus |f2T . 

Observing that [n]! for large n, is given by the g-Stirling approximation for q > 1 (see 
appendix for the explicit derivation) 



rr 



log [n]\ ~ n log [n] — — log q , (38) 

the entropy (|36|) can be written as 

c v- / i [TH + 9i] . i [ni +gi] . \ , Qn s 

S = <rii log + g i log — — rii g { \ogq\ . (39) 

This is similar to the structure of the entropy given by Eq.(^) and therefore the ex- 
tremization procedure can be carried out as was done before and derive the same condition 
as in Eq.(^) except for the factor g^. We observe, however, that the partition operation 
into subcells is not rigorously true in this context because of the nonextensive property 
(nonadditivity of the g-basic number) of the expression for the entropy in Eq.(^). For this 
reason, the mean occupation number derived from Eq.fl3"P|) is not rigorously proportional to 
the factor g^. The nonextensivity implies that the result for the mean occupation number 
is not entirely independent of the manner in which the energy levels of the particles are 
grouped into cells. 



VI. IDEAL Q-BOSE GAS AND Q-BOSON CONDENSATION 

We shall now proceed to study the thermodynamic behavior of an ideal g-Bose gas 
and the phenomenon of g-boson condensation. For a large volume (and a large number of 
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particles), the sum over all single particle energy states can be transformed to an integral 
over the energy, as follows 

9 V r°° 

=r- —j= / dxx^fix), (40) 

i V7T A d J0 

where x = (3e, e = p 2 /2m is the kinetic energy and A = h/ {2nmT) 1 ^ 2 is the thermal 
wavelength. 

We anticipate that the ground state will be associated with macroscopically large occu- 
pation number rather than a zero weight due to g-boson condensation. For this reason we 
need to isolate the ground state and include the contribution from all the other states in the 
integral. The number density of particles can thus be written as 

N 2 1 /-oo 1/2 1 , /z-V-l\ n 

77 = ~r 7? / dx x \ lo S — — + 77 > 41 

V V7T A 6 Jo iogq \z L e x — qJ V 

where no is the mean occupational number of the zero momentum state 

n = ^~ log f ^—^-) . (42) 
logg \l-qz) 

As in the standard boson case, we need to set the range of fugacity z which will correspond 
to non-negative occupation number. In the case of g-bosons we see that the condition is 
z<l/gforg>l and z < 1 for q < 1. It should be pointed out that we also have to 
require the existence of the JD of the mean occupation number which is encountered in the 
calculation of thermodynamic quantities such as the specific heat and this changes the upper 
bound of the fugacity z. We thus find the correct condition to be z < z q , where we have 
defined 

f q- 2 if q > 1 ; (A3) 
Zq ll if g < 1 . [ 6) 

We will have g-boson condensation when the critical combination of density and tem- 
perature occurs such that the fugacity will reach its maximum value z = z q . 

Following the prescription of the JD in the g-deformed thermodynamics derivatives, we 
obtain the expression for pressure above the critical point 



P 
T 

and below the critical point we have 



> = ^9 5/a (z,q), (44) 



P 
T 



< = 9 5/2 (z q ,q) ■ (45) 



Similar expression can be found for the number of particles above the critical point 



N 
V 



= ^3^3/2(^9). (46) 
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and below the critical point we have 

N 
V 



y + ^9 3/2 {z q ,q). (47) 



In the above equations we have defined the g-deformed g n (z, q) functions as 

( \ 1 n n-i 1 , (z- l e x -l\ 
11 logg \z L e x — q J 

5Z Un+l 5Z t,n+l ) ' 




logg V^i kn+1 t^i kn+ \ 

In the limit q — > 1, the deformed g n (z, q) functions reduce to the standard g n {z). In Fig. 
1 and 2 we present the behavior of g 3/2 (z,q) and g 5/2 (z,q) as a function of z for different 
values of the parameter q. 

The internal energy can be calculated considering the thermodynamic limit of Eq.(p6|) 
by means of the JD recipe. Using the expression for the pressure, Eq.(fP|), it is easy to 
verify that as in the undeformed case, the following well-known relation is satisfied for the 
g-bosons, 

U = | PV . (49) 

We can calculate the critical temperature by using the same method as in the standard 
boson case. Comparing the ratio of the critical temperature T c 9 of the g-deformed gas with 
that of the standard boson T c at the same density, we find 



TJ _( g, /3 (l) \ 2/3 
T c \9 3/2 (z q ,q)J 



(50) 



where <7 3 , 2 (1) = 2.61 is the value of the undeformed function when z — 1. In Fig. 3 we show 
the plot of the above ratio as a function q. We observe that the critical temperature of the 
g-boson is always higher than the standard boson and for q > 1 there is a rapid increase of 
the critical temperature T c 9 for small values of q. For example, for q = 1.01, T c 9 increases by 
18% and for q = 1.1, T c 9 increases by 75% with respect to the standard value. 

Applying the thermodynamic limit to the entropy of the g-boson in Eg. Q2D|) , we obtain 
the entropy per unit volume above the critical point with a structure similar to that of the 
standard boson, 



S_ 

V > = 

and below the critical point 



1/5 

^3 (2 95 ^ Z,q " > ~ 9 s/2 ( z ^) l °Sz) , (51) 



S_ 
V 



5 1 

2^3^6/2(^9) • (52) 



Let us observe that the generalized g-boson entropy obeys the third law of thermody- 
namics. In fact, in Eg. fl52"D , g B . 2 (z q ,q) has a finite value that depends on g and the entropy 
approaches zero in the limit of zero temperature. 
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As in the ordinary Bose condensation it is possible to show that in the g-boson conden- 
sation also a Clausius-Clapeyron equation holds and first order phase transition occurs. In 
fact it is easy to see that below the critical point the following equation is satisfied 



dP 
dT 



(53) 



where v c is the critical specific volume, defined as 

A 3 



(54) 



L q is the q-deformed latent heat given by 



TA ,_ 5 T 9 6/a (z q ,q) 
2 9 3/2 (z q ,q) 



(55) 



and where As is the difference in specific entropy across the transition region. 

We now proceed to calculate the heat capacity of the g-boson gas, starting from the 
thermodynamic definition 



_ dU 
v dT 



(56) 



V.N 



For this purpose we first need the derivative of the fugacity with respect to T (or ft), 
keeping V and N constant. To apply the JD prescription described earlier in Sec. IV, we 
start from the expression for the total number of particles, Eq. (|25|) and the identity (since 
the number of particles is kept constant) 



d 



dft^ log \T 



1 — Ki 



qK ij 







(57) 



where Ki = z e ^ . This identity can be rewritten according to our JD recipe as 

dKs r„\ , / 1 Ki 



Y 9ft K ' e 







l — qK i/ 

and now evaluating in the limit V — > oo, we obtain 

3 1 V^g B/2 (z,q) 



(58) 



idz 

~zdft 



V , N 2 ft V $g { 



(59) 



We shall now proceed to calculate the heat capacity. Using the discrete expression of 
internal energy (|26|), Eq.(|56|) can be expressed as 



c v = -ft 2 J2^ 



dKi 1 
dft logg 



log 



qKi 



(60) 
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Carrying out the limit V — > oo and utilizing Eqs.(^) and (0), we obtain the following 
expression for the specific heat per particle above the critical point 

15 zVtog r/a (z,q) 9 zV^g 5/2 (z,q) V^g 5/2 (z,q) 



~N 



(61) 



4 g 3/2 (z,q) 4 g 3/2 (z,q) vi q) g 3/2 (z,q) 
and similarly the specific heat below the critical point 

_ 15 zV^g 7/2 (z,q)\ z=Zq 
N < 4 X 3 /v 

where v is the specific volume below the critical temperature that can be expressed, by 
means of Eq . fl47D , in terms of the critical temperature T c as follows 

v 1 /T\ 3 / 2 



(62) 



A 3 9 3/2 (z,q)\Tr ■ (63) 

The above expressions have the same structure as that of the undeformed boson but the 
difference arises from the property zV^ q 'g n (z,q) ^ g n ^ 1 (z,q), where the equality is true for 
ordinary derivative only. From this observation it is easy to see that in the limit q — > 1 the 
specific heat reduces to the well-known undeformed result. 

As usual, the classical limit can be achieved considering the limit z — > 0. In this limit 
the deformed g n (z,q) functions reduce to 

o-l 

9n(z,q)^- z, (64) 

logg 

and from Eq . (|61~|) , the "classical" limit of the specific heat per particle number reduces to 

C v 
~N 



3q ~ 1 (65) 



d 2 logg 

As discussed before in the context of the entropy, the g-deformation persists also in the 
classical limit. 

We expect small deviations from undeformed behavior in the experimental observables, 
therefore only values of g close to standard value g = 1 are physically significant. Small 
deformation leads to a negligible departure from the high temperature limit of the specific 
heat but implies sharp deformation of the behavior of the specific heat in the range of the 
critical temperature. To exemplify this feature we plot in Fig. 4 the specific heat as a 
function of T/T c 9 for q = 1.05. We have chosen a value of q in the range q > 1 because this 
region appears particularly interesting with a higher critical temperature for small q (see 
Fig. 3). For this value of q the critical temperature T q is increased by 48% relative to the 
standard boson case. We observe that for q ^ 1, the specific heat shows a discontinuous 
A point behavior. This is a characteristic of g-deformation as has been observed in other 
investigations 

Using Eqs.flBTD and (|62l) we can calculate the jump in the specific heat A(C V /N) at the 
critical temperature as a function of q. In Fig. 5 we plot this behavior. We observe that the 
jump is an increasing function of q. 

Although the model that we have investigated is based on the Hamiltonian of noninter- 
acting particles, we note that the jump in the specific heat is of the order of the experimental 
value in the case of Bose condensation in 87 Rb atoms . 
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VII. CONCLUSION 



The outstanding problem in the theory of g-bosons has been the lack of a demonstration 
that the thermodynamic relations follow from the g-calculus framework. 

In this paper, we have shown that the whole structure of thermodynamics is preserved 
if the ordinary derivatives are replaced by the Jackson derivatives following the prescription 
described in Sec. IV. We establish a fully consistent set of relations between the thermo- 
dynamic functions (partition function, internal energy, mean occupation number) and this 
enables us to derive the entropy of g-bosons. The g-deformed entropy so obtained has been 
shown to follow from the deformed statistical weight as a consequence of the g-combinatorial 



calculus known in the literature pi] . This result represents a close connection between the 
quantum deformed algebra and the quantum statistical approach. 

The expression for the entropy is nonextensive because of the non-additive property (||) 
of the g-basic numbers. We find that for g 7^ 1, the entropy remains deformed in the classical 
limit as is also true of the other thermodynamic functions. This can be understood by ob- 
serving that the deformation arises from quantum groups but the nature of the deformation 
is inherently contained in the theory. A similar feature is found in the nonextensive Tsallis 



statistics and infinite statistics where the deformation persists in the classical limit ||T5l , [25 



In this framework, we have studied the basic properties of the ideal g-Bose gas in the 
thermodynamic limit and the phenomenon of g-boson condensation. We find that the critical 
temperature of the g-boson is always higher than that of the standard boson. The behavior 
of the specific heat exhibits a discontinuity at the transition point, which is in qualitative 
agreement, for values of q close to unity, with experimental data in the case of a dilute gas 
of Rubidium atoms [ 26 1 . 

We observe that for an ideal Bose gase the specific heat is continuous. On the basis of the 
Ginsburg-Landau theory of A-points, a discontinuous behavior of the specific heat implies a 
broken symmetry in the transition characterized by an order parameter. The deformation 
of the algebra in g-boson theory implies a broken permutation symmetry of the standard 



boson wave function. Therefore, the recent experimental data |26] can be interpreted as 
an indication of the effects due to g-deformation in Bose-Einstein condensation, where the 
order parameter of the phase transition depends on g. 

Although we employed the non-symmetric g-deformation in this investigation, all the 
results can be easily extended by using the symmetric g-calculus (g <-» g _1 ). We have 
confined our study to the g-deformation of bosons. It may be worthwhile to investigate the 
theory of g-fermions in this framework. 

Our theoretical framework and the results appear to provide a deeper insight into the 
behavior of the g-boson gas. We believe that the results derived here may be relevant to 
future investigations, and may be of interest from theoretical as well as experimental point 
of view. 
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APPENDIX A: 



Here we present a derivation for the approximation of g-basic factorial, [n]! for large n, 
which is the analog of the Stirling approximation. This is employed in the derivation of the 
entropy in Sec. V. We limit our discussion to the case of q > 1. 

Starting from the definition (f|) of the [x}\ and using the property (|T9| ) of the g-basic 
number, it is possible write any product factor contained in the [x]\ as follows: 



n 






[n\ 






n 


-1] 


= q~ 1 


[n] 






n 


-2] 


= q- 2 


[n] 






n 


-3] 


= q-* 


[n] 




- q- 2 - q^ 



[n-k] = q k [n] - q k - q k+l q 1 

[1] = q~ n+1 [n] - q- n+1 - q- n+2 q' 1 . (Al) 

For n>l, the leading term is seen to be 

[„]!« [n} n q-^V k (l-l) , (A2) 



where the first term arises from the product of the first term in each of the equations (ST) 
and the second term is the result of the sum of the dominant corrections. 
The above equation can be rewritten as 

[n]\ w [n] n q-^ 1 )' 2 ■ (A3) 

Taking the logarithm on both sides, we have 

9 

Tl 71 

log[n] ! w n log[n] - — logg - — , (A4) 
L q 

where we observe that the last term, which follows from the approximation: log(l — n/q n ) ~ 
—n/q n , is very small and significant only for q very close to unity (\q — 1| < 10~ 3 ). We 
neglect this term in the derivation of the entropy in Sec. V. 

We have verified numerically that the derived g-Stirling approximation is very good for 
large n. For example, for q = 1.5 it is correct to an error of 1.8% for n = 100, 0.2% for 
n = 1000 and 0.04% for n = 5000. 
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FIGURES 



gs/2 (z, q) 




FIG. 1. The behavior of g 3/2 (z, q) as a function of z for different values of q. The value q = 1 
corresponds to the standard g 3 , 2 (z) boson function. For q > 1 the upper bound of z is 1/q 2 and 
for q < 1 it is unity (see Eq.fl43|)) due to the existence of the JD of the g n (z, q). 



g5/2 (z, q) 




FIG. 2. Same as Fig. 1 for the function g 5/2 (z,q). 
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FIG. 3. The ratio T^/T c of the deformed critical temperature T% and the undeformed (q = 1) 
T c as a function of q. 



C v /N 




FIG. 4. The specific heat C v /N as a function of T/T« for q = 1.05. 
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q 



FIG. 5. The jump in the specific heat A(C V /N) at the critical temperature T c 9 as a function of q. 
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